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Distinct Scaling Regimes of Energy Release Dynamics in the 
Nighttime Magnetosphere 

1 1 2 .1 

V. M. Uritsky , E. Donovan , A. J. Klinias , and E. Spanswick 

I Based on a spatiotemporal analysis of POLAR UVI 

1 images, we show that the auroral emission events that 

initiate equatorward of the isotropic boundary (IB) ob- 

5^ ' 6 tained from a time-dependent empirical model, have sys- 

fvj . 7 tematically steeper power-law slopes of energy, power, 

, • 8 area and lifetime probability distributions compared to 

r—j ' 9 the events that initiate poleward of the IB. The low- 

H"~5 10 latitude group of events contains a distinct subpopulation 

pf-\ 11 of substorm-scale disturbances violating the power-law 

12 behavior, while the high latitude group is described by 

1 nearly perfect power-law statistics over the entire range 

I of scales studied. The results obtained indicate that the 

inner and outer portions of the plasma sheet are charac- 

^ ' 16 terized by substantially different scaling regimes of bursty 

O ' 17 energy dissipation suggestive of different physics in these 

i? , 18 regions. 
Oh. 

y i 1. Introduction 

C/2 ■ 

^^. 19 The activity of the nighttime auroral oval represents a 

(-H 20 wide range of dynamical processes in the magnetotail, 
C*) / 21 including substorm expansion onsets, pseudobreakups, 
I steady magnetospheric convection events with or with- 
1 out substorms, bursty bulk flows, and sawtooth events 
, (see e.g., Zesta et al. [2000]; Lui [2001]; Frey et al. [2004]; 
Henderson et al. \2QQfi\) . Despite the diversity of physical 
conditions associated with each particular type of auroral 
C^ I 27 activity, their net energy output can be described by a 
^^^ ' 28 set of apparently universal power-laws [Lui et al. [2000]; 
O ' 29 Lm [2002]; Untsky et al. [2003, 2002, 2006]) signaling 
I the existence of a organizing dynamical principle arrang- 
ing intermittent magnetospheric dissipation across vast 
I ranges of spatial and temporal scales. 
1 Power-law intermittency of energy dissipation has at- 

I tracted significant attention in modern statistical me- 
^ : 35 chanics (see Dhar [2006] and refs therein) and is often 
" "^ 36 considered a hallmark of turbulent and/or critical phe- 
/^ 37 nomena with no characteristic scales other than those 
V-j 38 dictated by the finite size of the system {Sreenivasan 
I et al. [2004]; Lubeck [2004]). Examples of such behav- 
I ior in geo- and space sciences include fully developed 
turbulence in hydrodynamic or magnetized flows {Lazar- 
I ian [2006]), Guttenberg-Richter statistics of earthquake 
1 magnitudes {Turcotte [1989]), scale-invariance in the so- 
I lar corona {Charbonneau et al. [2001]). In this context, 
the auroral activity provides one of the most impressive 
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46 examples of scale-free behavior in nature. The energy 

47 distribution of electron emission regions exhibits a power- 

48 law shape over a range of 6 orders of magnitude ( Uritsky 

49 et al. [2002]) which can be extended to up to 11 orders 

50 by combining the satellite data with ground-based TV 

51 observations {Kozelov et al. [2004]). 

52 The auroral emission statistics reported so far repre- 

53 sent global long-term properties of nighttime magneto- 

54 spheric disturbances. The fact that these properties are 

55 dominated by power-law scaling does not eliminate the 

56 possibility of a more complex behavior on the level of 

57 specific plasma sheet structures described by drastically 

58 different physical conditions and geometry. In this study, 

59 we are taking a step toward a better understanding of 

60 the relationship between the scale-free auroral precipita- 

61 tion statistics and the underlying central plasma sheet 

62 (CPS) morphology. We suggest that the inner and the 

63 outer CPS regions are responsible for three distinct scal- 

64 ing modes of the auroral precipitation dynamics, and pro- 

65 vide a possible physical interpretation for the observed 

66 differences. 

2. Data and Algorithm 

67 We have studied time series of digital images of night- 

68 time northern aurora (55-80 MLat, 2000 - 0400 MLT) 

69 taken by the Ultraviolet Imager (UVI) onboard the PO- 

70 LAR spacecraft in the 165.5 to 174.5 nm portion of the 

71 Lyman-Birge-Hopfield spectral band (integration time 

72 36.5 s, time resolution 184 s). The data analyzed in- 

73 elude 16,000 images covering two observation periods: 

74 01/01/1997 - 02/28/1997 and 01/01/1998 - 02/28/1998. 

75 Our analysis was based on spatiotemporal tracking of au- 
roral emission events {Uritsky et al. [2002, 2003]). The 
UV luminosity w{t, r) was studied as a function of time t 
and position r on the image plane. First, active auroral 
regions were identified by applying an activity threshold 
Wa representing a background UV flux. Contiguous spa- 
tial regions with w{r,i) > Wa were treated as pieces of 

82 evolving events. Second, by checking for overlap of com- 
es mon pixels between each pair of consecutive UVI frames, 

84 we constructed a set of spatiotemporal integration do- 

85 mains Ai(j = 1, .., A'') corresponding to each of A'^ individ- 

86 ual emission events found by our method. These domains 

87 of contiguous activity in space and time were used to com- 

88 pute the lifetime, T^, the energy, Ei = k J^ ■w{r,t) drdt, 

89 the peak power, Wt — fcmax(J' w{r,t) dr), as well 

90 the peak area. At — max( f dr) of every event, 

91 where k = 2.74 x 10^** J • photon^^ (Brittnacher et al. 

92 [1997]) and the integrals were numerically approximated 

93 by sums. The statistics reported below are for the thresh- 

94 old Wa — W photons • cm^ • s~^. Their main features re- 

95 main the same if the threshold is varied at least within 

96 the range 5 to 15 photons • cm^ • s^^. 

The auroral onset positions of each event were esti- 
mated with an error of about 300 km in either spatial 
direction. We organized the data to allow us to compare 
events that likely originate in the inner magnetosphere 
relative to events that likely originate in the outer mag- 
netosphere. To accomplish this, we determine the event 
locations relative to the isotropic boundary (IB) at the 
meridian of the event origin ^ : 

(j)i = MLati - [Ao - Al cos {jviMLTi ~ MLTo)/12)] (1) 

97 Here, MLTi and MLati are the onset coordinates of the 
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j-th event, ^o, Ai, MLTq are the coefficients of the em- 
pirical IB model due to Gvozdevsky and Sergeev [1995] 
which were computed based on hourly values of auroral 
electrojet index (WDC for Geomagnetism, Kyoto) and 
solar wind dynamic pressure (ACE spacecraft) at the be- 
ginning of each auroral intensification. The ions main- 
tain their isotropy on the poleward side of IB due to the 
effective pitch angle scattering in the tail current sheet 
{Newell et al. [1996]). As a result, the latitude of IB has 
a high correlation with the magnetic field inclination at 
the geomagnetic equator as well as with the equatorward 
boundary of the proton aurora as determined from di- 
rect optical observations {Donovan et al. [2003]). Here, 
we have used the empirical IB model to separate inner 
and outer emission events. Positive values of (j) on av- 
erage indicate events originating on field lines that map 
further out on the plasma sheet, as opposed to negative 
values corresponding to the events that originate closer 
to Earth. 

The statistics of the emission events that initiated on 
the poleward {(f) > 0) and on the equatorward {(j) < 0) 
side of the IB are characterized by sets of probability 
density distributions p{x), where x £ {E,W,A,T}. The 
power-law exponents obtained from these distributions 
are denoted as Tx, with the subscript indicating the vari- 
able under study. 

3. Results and Discussion 

In the end we have 7481 events, 6231 and 1250 pole- 
ward and equatorward of the IB, respectively. These 
yielded a database of parameters E, W , A, T, MLT, 
MLat, (j) for those events. Figure 1 is the emission en- 
ergy _E as a function of the relative onset latitude (j) for 
all events (n = 7481). The scatterplot has clear asym- 
metric shape suggesting different statistics for negative 
and positive (j> values. As we discuss below, these statis- 
tical subsets are characterized by significantly different 
regimes of scaling behavior indicating different physical 
environments. 

Figure 2 shows probability distributions for the events 
which initiated on the poleward and equatorward sides 
of the IB. The high-latitude (HL) events (0 > 0) ex- 
hibit stable power-law distributions of emission energy 
E, peak emission power W , lifetime T and peak area A 
(not shown) with near-zero regression errors. The ranges 
of scales of these power-law behaviors involves both small 
auroral activations and large events whose energy output 
lies in the range of global substorms. The exponents te, 
Tw, TA and TT are close to the values obtained earlier for 
the same observation period without filtering the activity 
by the onset location ( Uritsky et al. [2002] ) . 

The low-latitude (LL) events {4> < 0) have more com- 
plicated statistical properties as can be expected from 
Fig. 1. Their distribution functions demonstrate a 
crossover behavior involving small-scale regions described 
by Tx exponents which are greater than the correspond- 
ing exponents of HL events, as well as large-scale regions 
where the slopes are significantly shallower (the Tx expo- 
nents for these regions are not well defined due to the in- 
sufficient number of large-scale LL events) . The position 
of the crossover E* = 5 x 10^^ J on the p{E) distribution 
of LL events matches the data gap in Fig. 1 separating 
low- and high-energy events with negative relative lati- 
tudes (j). 

As can be seen from Table 1, only 1/6 of the events 
contribute to the population with < 0. However, their 
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participation in the energy budget of the electron aurora 
approaches 50 percent as can be seen from the compar- 
ison of the overall energy Etot released by LL and HL 
events. The fact that the HL events are the the most 
common form of emission dynamics (relative occurrence 
frequency 82.3%) is not a surprise as these events likely 
map to the reconnection regions of the plasma sheet, as 
we briefly discuss below. The table also shows that the 
scaling behaviors of LL and LH events are described by 
subsantially different sets of r exponents. 

Table 2 provides scaling exponents of energy distri- 
butions p{E) as well as several parameters reflecting 
the state of the solar wind and the auroral magneto- 
sphere at the beginnings of emission events for four non- 
overlapping ranges of magnetic latitude MLat. The data 
show that the group of events that initiated southward of 
MLat=65 are characterized by significantly higher statis- 
tical values of solar wind dynamic pressure and electrojet 
index. These events seem to appear in a rather stretched 
magnetotail configuration as refiected by the magneto- 
tail (MT) index {Gvozdevsky and Sergeev [1995]). One 
can also notice that the lower the magnetic latitude, the 
higher is the exponent te- This tendency is in accordance 
with the difference in the LL- and HL-event statistics re- 
ported in Table 1, and it also suggests that the change 
of the scaling regime of auroral emission events with on- 
set latitude is a continuous process rather than a discrete 
transition fully dictated by the empirical IB location. 

From the statistical mechanics viewpoint, the exis- 
tence of two or more scaling regimes within the same 
physical system signals the presence of several distinct 
universality classes (UC) in the underlying turbulent 
dynamics (Luheck [2004]; Dhar [2006]). It is known 
that many stochastic scale-invariant phenomena can be 
mapped onto a finite collection of such classes. This map- 
ping is insensitive to a variety of physical parameters and 
conditions and can be the same for quite different sys- 
tems. Typically, each UC brings its own set of critical 
exponents which express its inherent set of relevant vari- 
ables and the underlying small- and large-scale symme- 
tries {BenHur and Biham [1996]). In this context, our 
results clearly show that the energy release dynamics in 
the inner and outer CPS belong to different UCs and are 
therefore governed by different physical mechanisms. 

The nature of the isotropic boundary sheds some light 
on the scaling regimes observed. On average, HL events 
tend to initiate in the outer CPS region. The energy con- 
version in this region is believed to be dominated by mag- 
netic reconnection {Birn and Hones [1981]). Our results 
strongly suggest that near-Earth (midtail) reconnection 
is a turbulent bursty process with no well-defined dissi- 
pation scales. Judging by the relative occurrence of these 
events (Table 1), the scale-free activity in the reconnec- 
tion region is the dominant mode of energy release dy- 
namics in the nightside magnetosphere. The power-law 
emission exponents obtained for this region are rather 
close to the exponents from a driven current sheet simula- 
tion (Klinias et al. [2004]) consistent with multiscale tur- 
bulent reconnection being the source of scale-invariance 
in the outer CPS dynamics. The LL events are mainly 
produced in the inner CPS. With a less stretched mag- 
netic field topology, this region is not a region that favors 
magnetic reconnection; however, it can be prone to cur- 
rent disruption which offers an alternate mechanism for 
energy release in the inner tail [Lui [2000]). 

Our findings can be summarized as follows: (1) the 
high-latitude emission events are characterized by broad- 
band power-law statistics with no characteristic scales; 
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230 (2) the low-latitude events constitute a non-uniform sta- 

231 tistical population with an energy crossover separating 

232 large and small scale activity; (3) the scaling exponents 

233 of the low- and high-latitude events are significantly dif- 

234 ferent within the entire range of scales studied, which 

235 signals distinct universality classes of energy release dy- 

236 namics. To the extent that the empirical IB model can 

237 be considered an accurate representation of the true IB, 

238 the onsets of the HL events are located in an outer CPS 

239 region, the onsets of the LL events in an inner CPS re- 

240 gions, and these two regions exhibit distinctly different 

241 energy release dynamics. 

242 The difi'erence in UC following from this analysis is 

243 suggestive of different physics. It would be a significant 

244 step forward if we could use physics-based models of CPS 

245 dynamics to predict UCs, and to thereby determine how 

246 the different statistical behaviors of the inner and outer 

247 regions might be generated. 
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Notes 

1. A more accurate analysis based on Tsyganenko T96 and 
T05 field-line models as well as optical determination of 
the isotropic boundary will be published elsewhere. 
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Figure 1. Scatterplot of emission energies E versus 
relative magnetic latitudes (j> of onset locations mea- 
sured with respect to the empirical IB latitude (IBL) 
(see eq.(l)). The high-latitude {<j} > 0) and the low- 
latitude {<j} < 0) events contribute to two distinct scaling 
regimes of energy release dynamics in the nighttime mag- 
netosphere as discussed further in the text. The dashed 
horizontal line marks the position of the crossover on the 
energy distribution shown in Fig. 2. The sharp lower 
cutoff is due to the fixed activity threshold Wa- 
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Figure 2. Probability distributions of energy E, peak 
power W , and lifetime T of auroral emission events initi- 
ated poleward {(j) > 0, crosses) and equatorward {4> < 0, 
stars) of the empirical IB. The low-latitude distributions 
are shifted for easier comparison. The solid lines show 
the log-log distribution slopes for small-scale LL events 
(the power-law portions to the left of the dashed vertical 
lines) and for the entire group of HL events. 



Table 1. Comparative cliaracteristics of low- and iiigh- 
latitude emission events. 



Type of events 


LL (</. < 0) 


HL (0 > 0) 


# of events 
Etot,J 

Te 

TW 
TT 
TA 


1250 (16.7%) 
5.7 X IQis (47.5%) 

1.96 ±0.04 
2.14 ±0.07 
3.30 ±0.26 
2.30 ±0.10 


6231 (82.3%) 
6.3 X 1015 (52.5%) 

1.56 ±0.02 
1.80 ± 0.02 
2.26 ±0.13 
1.93 ±0.06 
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Table 2. Scaling exponent te of emission events, solar wind 
dynamic pressure {Pj, auroral electrojet (AE) index and the 
magnetotail (MT) index (IB latitude at MLT=0, (eq. 1)), for 
four ranges of onset MLat. 



MLat 


TE 


Pd, nPa 


AE, nT 


MT 


<65 
65.. .67 
67.. .69 
>69 


2.06 ± 0.04 
1.88 ±0.04 
1.59 ±0.03 
1.58 ±0.03 


2.64 ± 0.05 
2.42 ± 0.03 
2.41 ± 0.03 
2.20 ±0.03 


206 ±5 
106 ±3 

84 ±3 
93 ±4 


63.2 ± 0.03 
63.9 ± 0.03 

64.0 ±0.03 

64.1 ± 0.03 



